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Optical properties of large arrays of isolated quantum dots 
are discussed in order to interpret the existent photolumines- 
cence data. The presented theory explains the large observed 
shift between the lowest emission and absorption energies as 
the average distance between the ground and first excited 
states of the dots. The lineshape of the spectra is calculated 
for the case when the fluctuations of the energy levels in quan- 
tum dots are due to the alloy composition fluctuations. The 
calculated lineshape is in good agreement with the experi- 
mental data. The influence of fluctuations of the shape of 
quantum dots on the photoluminescence spectra is also dis- 
cussed. 



I. INTRODUCTION 

Reduced-dimensionality .-structures are currently at- 
tracting much attentionJjta Modern technology made it 
possible to create nanoscale OD structures where motion 
in all directions is quantized (quantum dots). Often a 
system represents a large array of independent quantum 
dots. Because of extremely small dimensions, the fluctu- 
ations of parameters of individual quantum dots become 
an important factor since they determine the properties 
of the whole array. 

This paper was initiated by the experimental work Ref. 
[j], where nice photoluminescence (PL) and photolumi- 
nescence excitation (PLE) experiments were performed 
with an array of self- assembling quantum dots. The con- 
ventional nonselective (with above-barrier excitation) PL 
reveals a broad peak of about ~ 50 meV in width. This 
width originates from a wide spread of energies of differ- 
ent quantum dots in the array. Both PLE and selectively- 
excited PL spectra (when only the quantum dots that are 
in resonance with incident light are excited) show sets of 
broadened peaks with 2-3 times smaller widths. These 
peaks correspond to the distribution of energy levels in 
the subset of quantum dots that are resonantly excited. 
It was observed that the first PLE peak is strongly shifted 
from the detection energy. The origin of this shift re- 
mained unclear. Moreover, no measurable Stokes shift 
was observed in a recent paper Ref. ||, where the PL has 
been measured together with the direct absorption by the 
layer of quantum dots. 

In this paper a simple theory of the PL from an ar- 
ray of quantum dots is developed. We suggest that 
the process of photoexcitation of a dot into its lowest 
optically-excited state does not contribute to the PL 



signal. The proposed interpretation of the experimen- 
tal data explains large Stokes-like shift between the PL 
and PLE peaks as the average distance between the two 
lowest optically-excited states. In what follows we use 
the term "ground state" to denote the lowest optically- 
excited state because it corresponds to the ground state 
of the photoexcited carriers in the dot. 

It is suggested that the PL and PLE lineshapes are 
completely determined by the statistical distribution of 
the energy levels of different quantum dots in the array. 
More precisely, it is determined by the distribution of 
pairs of energy levels. It is shown that such distribution 
is essentially two-variable. That is, there exists a corre- 
lation in the positions of different energy levels in each 
quantum dot throughout the array, however, such corre- 
lation is not 100%. We show that this feature causes the 
difference in the positions of the maxima of the PL and 
PLE spectra. 

The fluctuations of energy levels due to the random 
potential caused by the alloy composition fluctuations 
are studied in Section IV. We show that the major part 
of the observed linewidth can be accounted for by this 
mechanism. We also suggest that the first two excited 
states observed in Ref. [0 originate from the twofold de- 
generate first excited level, when the degeneracy is lifted 
by the random potential. The density of states and the 
two-level distribution function, which accounts for the 
correlation in energy-level positions, are calculated. The 
spectra determined by these functions describe most of 
the features of the experimental spectra. 

Finally, the effects caused by fluctuations of the shape 
of the quantum dots are discussed. It is suggested that 
the fluctuations of the shape of the dots should increase 
the energy-level correlation. 



II. ORIGIN OF THE STOKES-LIKE SHIFT 

In the present paper we consider the case when both 
the electron and the hole are confined in the quantum 
dot. Energy of the quantum dot is measured from the 
unexcited state of the dot with no electrons and holes. 
The term "ground state" refers to the ground state of 
the quantum dot with an excited electron-hole pair, i.e. 
to the lowest optically-excited state. 

The key feature of an array of quantum dots that dif- 
fers it from a quantum well or bulk material is that there 
is no charge transfer between the dots, or at least such 
transfer is strongly suppressed. Different quantum dots 
therefore give independent additive contributions in any 
optical experiment .□ 
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The density of states for each quantum dot represents a 
set of delta-function-like peaks, while the total density of 
states of the whole system can be spread in a wide energy 
range due to the inhomogeneous broadening. Several ex- 
periments have been reported recently where the con- 
tributions to the luminescence from single quantum dots 
were found.cTQQ Single quantum dots give extremely nar- 
row sub-meV spectral lines. When, however, the number 
of the excited quanturLL-dots is large, a broad lumines- 
cence peak is observed. Btrtl 

We suggest that such a feature makes it difficult to 
probe optically the ground states of the dots and thus 
causes an apparent Stokes-like shift observed in Ref. [I]. 
Probing of the ground states of the dots may require 
special technique like high resolution or time- resolved PL. 

Indeed, let's consider the contribution to the PL from 
the process when the photon is absorbed into the ground 
state of a dot and then re-emitted. If the optical process 
is not phonon- assisted, the energies of the incoming and 
outgoing photons are exactly equal. The PL response 
can hardly be observed since it is hidden by the incident 
beam scattered by other elements of the experimental 
environment. This situation is quite different from that 
in quantum wells or bulk semiconductors, where an elec- 
tron excited to the conduction band has always quantum 
states with smaller energies. It can loose some energy 
before it recombines in an independent optical process. 

In a quantum dot, the emitted photon may have its en- 
ergy below (or above) that of the exciting light if the PL is 
phonon-assisted. However, the phonon must be emitted 
(absorbed) together with photon in the same quantum 
process. The probability of such process is determined in 
higher-order perturbation theory in the electron-phonon 
coupling constant and is therefore much smaller than the 
probability of the direct transition. 

In order to give a substantial contribution to the PL 
signal, the dot must be pumped into one of its excited 
states. Therefore, the minimum distance between the 
excitation and detection energies seen in the spectra is 
equal to the distance between the ground and first ex- 
cited states of the dot. In fact, no Stokes shift was found 
between the PL and the direct absorption by a layer of 
quantum dots measured in Ref. [|. 

III. POSITIONS OF PEAKS IN THE PL AND 
PLE SPECTRA 

In this section the shape of the PL and PLE spectra 
is described qualitatively. Different peaks observed in 
the spectra are assigned. It is shown that existence of 
a random spread of interlevel distances in quantum dots 
causes substantial deviations of the positions of peaks 
seen in PLE and selectively-excited PL spectra. 

We shall assume for simplicity that the energy relax- 
ation in a quantum dot occurs faster than recombination, 
so that the light is always emitted from the ground state 
of the dot. 



Let us first consider a simple model when the distances 
between different energy levels are the same for all quan- 
tum dots, however, there is a wide distribution of energy 
levels in the array. In other words, let the picture of en- 
ergy levels be the same in each quantum dot but shifted 
randomly as a whole. This implies a 100% correlation in 
the positions of different energy levels in each dot. 

If the conventional PL technique is used so that the 
pumping is performed with energies well above the bar- 
riers between the dots, all the dots are excited and emit 
light at their ground-state energies. The emission then 
represents a broad peak due to the wide distribution of 
the ground-state energies. 

Let us now consider PLE and selectively-excited PL 
measurements. Here only those dots are excited, which 
have one of their energy levels in resonance with the excit- 
ing light. As suggested above, the dots pumped into their 
ground states do not contribute to the spectra. When the 
interlevel distances are the same in all dots, both PLE 
and selectively-excited PL spectra would show a set of 
delta-function- like peaks as sketched in the Fig. 1. 

It is more convenient to start with the PLE (left). By 
fixing the detection energy one selects the subset of all 
quantum dots in the array with the ground-state energy 
Eq = Edetector (we have assumed that the light is always 
emitted from the ground state of the dots). The PLE 
signal appears when the excitation energy matches the 
energy of an excited state {Ei aser = E\, E%, etc.) in 
the selected subset. Thus the first observed (lowest in 
energy) PLE peak corresponds to the first excited energy 
levels of such dots that have their ground-state level at 
the detection energy. 

The analysis of the selectively-excited PL (right) is 
somewhat more complicated but similar. The fixed en- 
ergy of excitation selects several subsets of all quantum 
dots such that E\ = Ei aser , E 2 = Ei aser , etc. A nonzero 
PL signal appears when the detection energy matches the 
ground-state energy of one of the subsets. The position of 
the first observed PL peak (highest in energy) thus corre- 
sponds to the ground state energy of such dots that have 
their first excited energy level at the excitation energy. 
The distance between the first two PL peaks is equal to 
E21, the distance between the first and the second excited 
states of the dots. 

When interlevel distances are the same in all dots, all 
distances between different peaks in the PL (PLE) series 
also remain the same, while the whole picture shifts with 
the shift of the laser (detector) energies. It means that 
positions of these peaks being plotted against the laser 
(detector) energies must form a set of straight lines with 
the unit slope. 

The slopes of such lines obtained from the PLE and 
selectively-excited PL spectra obtained in Ref. [l] are all 
less then 1. For the first PL and PLE peaks, rf-g-, 

dEmax.l/dElaser = 0.91, dE^^/dE detector = 0.77-lil 

It is easy to see that such deviation cannot be at- 
tributed to the dependence of the interlevel distances on 
the ground-state energy. Indeed, in this case the slopes 
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obtained from the PL and PLE spectra should be inverse 
of each other. We see, however, that 

clF PL ( dF PLE \ _1 

max.m I max.m j 

dEi aser \ dEdetector J 

We show below that existence of a spread of interlevel 
distances in the array causes deviation of the positions 
of peaks seen in the PL and PLE spectra in such a way 
that both derivatives in Eq. (1) become less then 1. 

To demonstrate this, let us assume that the interlevel 
distances in the dots fluctuate randomly in some (nar- 
row) energy interval. In this case selection of a subset of 
quantum dots by fixing one of the energy levels does not 
yet determine the positions of other energy levels. Then 
one should observe a sequence of broadened peaks in both 
PL and PLE. The scale of broadening is determined by 
the distribution of interlevel distances only, therefore it 
can be narrower than the distribution of the absolute en- 
ergies. As shown below such behavior is natural if the 
spread is caused by a random potential. 

To understand how the positions of maxima of the 
broadened peaks are shifted, it is convenient to draw a 
three-dimensional picture shown in Fig. 2. Here the in- 
tensity measured by the detector is plotted as a function 
of two variables — the excitation and the emission ener- 
gies. Clearly, such a plot contains all information, which 
both PL and PLE can provide. To get the shapes of 
the PL or PLE spectra one simply has to slice the three- 
dimensional plot along the Ei aser or Edetector axis. 

First, the intensity is zero in the half-plane Ei aser < 
Ede.te.ctor- The intensity is nonzero only when Ei aser and 
Edetector are in resonance with the excited and ground 
states of some dot simultaneously. If the spread of inter- 
level distances is small, the three-dimensional plot has a 
shape of a set of narrow ridges, elongated in the direction 
parallel to the line Ei aser = Edetector- Each ridge corre- 
sponds to the optical process for which Edetector = Eq 
and Ei aS e r = E m , m=l, 2, etc. The width of each ridge 
is determined by the spread of the corresponding inter- 
level distances, while the length is larger and represents 
the large spread of the ground-state energies. The dis- 
tance between the ridge and the line Ei aser — Edetector is 
determined by the average distance between the ground 
and the corresponding excited states of the dots. 

The inset in Fig. 2 shows a fragment of the same plot in 
the isoline projection (dashed lines). Crosses and pluses 
show the positions of the maxima in the PL and PLE 
crossections respectively. The "PL and PLE lines" are 
the loci of the maxima observed in the PL and PLE spec- 
tra. 

As one can see, the positions of the maxima deviate 
from the major axis of the ridge and from each other. In- 
deed, these positions lie in such points where the crossec- 
tion line (parallel to one of the axes) is tangential to the 
lines of equal intensity. Note that the deviation of the 
positions of the maxima is such that both derivatives in 
the Eq. (1) are less than 1. 



Both lines intersect exactly at the maximum of the 
ridge, which coordinates give the average positions of the 
ground and the excited states of all the dots in the array. 

In general, if the oscillator strengths are the same 
for all allowed transitions, the intensity as a func- 
tion of the excitation and detection energies is propor- 
tional to the mutual level-level distribution function: 

I (Edetector : E[ aS er^) ^2/m ^ (.-^detector Eo,E[ ase r — 

E m ), where E m are the energy-level positions, E be- 
ing the position of the ground state. In the following 
section we derive the shape of the distribution function 
P(Eq, E m ) when the spread of the energy- level positions 
is caused by the composition fluctuations. We show that 
a quite complicated PL and PLE lineshape observed in 
Ref. |l| can be sufficiently well described in this way. 

IV. FLUCTUATIONS OF ENERGY LEVELS IN 
QUANTUM DOTS 

In this section we study the properties of the statistical 
distribution of the energy levels in quantum dots caused 
by a white-noise random potential. We present a strong 
evidence that the major part of the observed spread of 
the PL and PLE peaks is caused by universal composi- 
tion fluctuations in the dots. These fluctuations are a 
generic property of semiconductor alloys and produce a 
theoretical limit for unification of quantum dots in an 
array. 

We also suggest that the first two excited levels that re- 
veal themselves in the PLE and selectively-excited PL ex- 
periments originate from a single doubly degenerate level 
when the degeneracy is lifted by the random potential. 

In semiconductor alloys the lattice sites are occupied 
randomly with two types of substitutional atoms. We ig- 
nore here the correlation between occupation of different 
sites. The composition x averaged over a small volume 
always fluctuates. The order of magnitude of the fluctu- 
ations is inversely proportional to the square root of the 
volume over which the averaging is performed. Though 
usually small, this effect can be important in extremely 
small structures. 

In small quantum dots such composition fluctuations 
cause shifts of energy levels from the positions deter- 
mined by the average composition xq. The "local" com- 
position x varies from dot to dot, and also inside the dot. 
According to that, the energy levels fluctuate from dot to 
dot. Shifts of the energy levels in different quantum dots 
are independent from each other. Inside a single quan- 
tum dot the shifts of different energy levels are correlated. 
Note, however, that such correlation is not 100% (as it 
would be if there was only one fluctuating parameter like 
the diameter of a dot). 

In order to calculate the statistical distribution of en- 
ergy levels caused by composition fluctuations we use the 
method developed by Efros and Raikh (for a review see 
Ref. uu . This method is applicable if the size of the wave 
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function is much larger than the lattice constant. In this 
case one needs to know only the shape of the wave func- 
tion in a quantum dot and the slope of the dependence of 
the gap on composition, dE g /dx, at average composition 
x . 

The result depends on whether the unperturbed en- 
ergy levels are degenerate or not. Without degeneracy 
the distribution of energy levels E m is Gaussian with the 
standard deviation given by 



7 / d\ V4(r), 



(2) 



where e m = E rn — E m , ip m is the wave function (real) 
corresponding to the energy level E m , and 7 is given by 
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dE g \ x {l-x ) 



dx 



N 



(3) 



where N is the number of lattice sites per unit volume. 

The covariance between e m and e n can be obtained in 
a similar way: 

e^T = pa m a n = 7 J d 3 r ^(r)^(r). (4) 

The coefficient p, p < 1, is called the coefficient of corre- 
lation between E m and E n . 

To find the shape of the PLE and selectively-excited 
PL spectra we need the mutual distribution function for 
Eq and E m . The most general form of the two- variable 
Gaussian distribution is given by 



G 2 (e , ;C0,cr m ,/j) 



1 



2ir<j <j m y/l - p 2 



x cxp 



r ^2 



2(1 -P 2 ) 



a aocr„ 



(5) 



This is just an analytical expression for the shape of the 
ridge, discussed in the previous section. The ridge is 
strongly elongated when the parameter p is close to 1. It 
is equal to 1 in the limiting case when an exact relation 
between the energy-level positions exists, so that the two- 
variable statistical distribution becomes effectively one- 
variable. The ratio a m /<Jo determines the orientation 
of the ridge. The ridge is parallel to the line Ei aser = 

Edetector when <J m = (To. 

The situation is, however, more complicated if a degen- 
eracy exists. The random potential shifts the degenerate 
energy level and lifts the degeneracy. The distribution 
of energies in the vicinity of an unperturbed degener- 
ate level appears to be not Gaussian. The mutual two- 
variable distribution function for a transition between E 
and E m is not Gaussian either. Instead, it has a shape 
of two close parallel ridges corresponding to each of the 
split-off energy levels. 

For simplicity we restrict ourselves to the axially sym- 
metric quantum dots, where each energy level except the 



ground state is doubly degenerate. In this case it appears 
to be possible to derive the general form of the distribu- 
tion function without knowledge of the shape of the wave 
functions in the quantum dot. 

For an axially symmetric system two degenerate wave 
functions with the angular momentum |m| have the form 
ip m ±(r) — ip m {r)e ±lm ^ , where tjj m can be made real. The 
positions of energy levels E m ±, split and shifted by ran- 
dom potential of the particular configuration, can be ob- 
tained as the eigenvalues of the secular matrix 



SH = 



u x + iy 
x — iy u 



(6) 



where the matrix elements u and x + iy take random 
values in each quantum dot and are given by: 



u = yfi j d 3 r V(r)V4(r) 
x + iy = V7 J d 3 r V(r)^(r)e 



2im<p 



(7a) 
(7b) 



Here V(r) is the white-noise random potential with cor- 
relator < V(r)V(r') >= 7(5(r - r'). 

Eigenvalues of SH are e± = u ± yjx 2 + y 2 (the energy 
is measured from the unperturbed energy level). It is 
easy to see that u, x, and y are independent Gaussian 
random variables with u 2 — cr^, x 2 = y 2 = er^/2, o 2 m 
being given by Eq. (2). 

The density of states in the vicinity of E m (the distri- 
bution function for e m = E m — E m ) is given by 

P{tm) = ^2 J J J dudxd V S(e m -uT V x2 + V 2 ) 



xGi(u; cr m )Gi(a;; ^=)Gi(y; ^=) 
= D (e m ; a m ), 



(8) 



where G\ is the standard one-variable Gaussian distri- 
bution, die- a) = a^die/a-l), Gi(z;l) = (2tt)- 1 / 2 x 
cxp(— z 2 /2); and we have introduced the notation for the 
distribution function D (e;<r) = a~ 1 D (e/<T; 1), 



D 0{z ;l) = \f-e- 



2 /2 + 



3V3 



erf(^). (9) 



The meaning of the index in D (e; a) will soon become 
clear. 

The Eq. (9) describes a bell-shaped curve, which de- 
termines the density of states and, hence, the absorp- 
tion spectrum in the vicinity of the excited level E m . It's 
shape can be fairly well approximated by a Gaussian with 
the effective dispersion a e ff = a m ^/2, as shown in Fig. 
3. 

In order to find the mutual two-level distribution func- 
tion that determines the shape of the PL and PLE spec- 
tra it is important to account for the correlation between 
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matrix elements of SH and the shift of the ground-state 

It is useful to note that the expression for the ma- 
trix clement u is the same as the expression for the shift 
of a non-degenerate level with the wave function ipmij')- 
The value of u is, therefore, correlated with eo via the 
same two- variable Gaussian distribution as in Eq. (5): 
P(eo,u) = G2(eo,u; ao, a m , p), with the parameters uq, 
a mi and p given by Eqs. (2), (4). The parameters x and 
y are statistically independent of eo and u, because the 
corresponding covariances become zero when the integra- 
tion over the azimuthal angle (f> is performed. Thus, we 
obtain: 

-P(eo, £m) =^2 J J J dudxdy S(e m - u T \J x 2 + y 2 ) 
xG 2 (e ,u;a ,a m ,p)G-i(x; ^)G x {y; 



= Gi(e ; &o)Dp(e m — e p— ; cr m ), 
where the function D p (e;a) = a -1 D p (e/o~; 1), and 
D^;l) = ^^exp / 



(10) 



y/wp 

2z 



p-1 



3/2 



exp 



VW-1) 



(11) 



Here p = 3 — 2p 2 . The function Dq(c; a) defined by Eq. 
(9) is a particular case of (11) with p = 0. 

The function D p (e m ] a m ) gives the distribution of the 
energy sublevels in the vicinity of E m , when the position 
of the ground state is fixed. It is normalized in such a 
way that J de D p (e;a) = 2 (according to the twofold 
degeneracy), and J de e 2 D p (e;a) — 2a 2 (2 — p 2 ). The 
function D p (e; a) is symmetric in e. It has one maximum 
when p < l/v2 and two maxima, when p > l/y/2. The 
maxima become more pronounced when p tends to 1. 

The function P(e ,e m ) determined by Eq. (10) gives 
the probability density for a quantum dot to have the 
ground state at the energy Eq + eo and an excited state 
at the energy E m + e m . It is proportional to the intensity 
measured by the detector at fixed excitation and detec- 
tion energies. Hence, it determines the shape of both PL 
and PLE spectra when the proper argument is fixed. 

The function P(eo, e m ) appears to be not very sensitive 
to the ratio o~ m /<Jo- It is, however, quite sensitive to the 
value of the correlation coefficient p. This function is 
plotted in Fig. 4 for a m — o~q and two different values 
of p. As shown in the next section, it is natural for the 
coefficient p to be close to 1. The value p = 0.94 [Fig. 
4(a)] gives the best fit to the experimental data. When p 
tends to 1 [Fig. 4(b)] the intensity in the dip between two 
maxima approaches zero. The lineshapes of the spectra 
described by Eq. (10) are discussed in detail in the next 
section. 



V. DISCUSSION 

We replot the positions of the PL and PLE maxima 
observed in Ref. [I] on the combined plot in Fig. 5 to illus- 
trate that the experimental behavior is in agreement with 
our consideration (compare with the inset in Fig. 2). The 
dotted line shows unit slope, Ei aser = Edetector- Three 
lines correspond to the positions of two maxima observed 
in each PL curve (crosses) and one — in each PLE curve 
(pluses). Positions of only one PLE maximum for each 
Edetector are shown because the second PLE maximum is 
seen not clear enough. The constant value of 1290 meV 
has been subtracted from all energies. 

Three lines shown in Fig. 5 are consistent with the 
qualitative picture of two close parallel ridges. As sug- 
gested above, these ridges correspond to the two optical 
processes where E detector = Eq and Ei aser = E\±; E 
being the ground-state energy and E\± — the energies 
of the first and second excited states. 

We use the notation E\± instead of E\ t 2 for the two ex- 
cited states according to the idea that they originate from 
the twofold-degenerate first excited energy level E\ when 
the degeneracy is lifted by random potential. This idea 
is supported by the fact that the ratio of interlevel dis- 
tances appears to be (Ei + — Ei^)/(Ei- — Eo) ~ 0.6. For 
a cylindrical quantum well with infinite potential walls 
the corresponding ratio is about (E 2 i/E w ) = 1.3. 

The point of intersection of two lines gives the posi- 
tion of the maximum of the first ridge. The difference 
of 55 meV in its coordinates is nothing but the observed 
Stokes-like shift between emission and absorption ener- 
gies. The distance between two ridges measured along 
the Ei aser axis gives E\+ — E±-, the mean splitting of 
the excited state E±. 

Using Fig. 5 we may conclude that the average distance 
between the ground and (split) first excited state is about 
75 meV, while the average splitting of the excited state 
is approximately 40 meV. 

Fig. 6 shows the comparison between the presented 
theory and the experiment. The experimental spectra 
from Ref. [l] are replotted in Fig. 6(a). The theoretical 
PLE and selectively-excited PL spectra [Fig. 6(b)] are 
obtained with the use of Eq. (10) by fixing eo or e m re- 
spectively. Experimental values for the average positions 
of the energy levels Eq = 1276 meV and E\ = 1351 meV 
are used. 

The Fig. 6(b) is essentially the Fig. 4(a), replotted in 
the wavelength scale for better comparison. The value 
p = 0.94 gives the best fit for the data from Ref. [I]. It is 
seen that the curves obtained reproduce all the features 
of both PL and PLE spectra. 

The structure of the expression (10) is such that if eo 
is fixed, it describes a curve symmetric around the point 
e m = eopa m /ao. Thus (if p > l/\/2) the PLE line should 
reveal two symmetric peaks at any detection energy. This 
is exactly the behavior seen in the experiment. When the 
detection energy is changed, the whole spectrum must 
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shift linearly with it. Indeed, Fig. 5 shows that the po- 
sitions of PLE maxima depend linearly on the detection 
energy. The slope of the PLE curve in Fig. 5 gives the 
ratio cti/ctoj which appears to be close to 1. 

If the position of the excited level, e m , is fixed, the 
lincshape is asymmetric. When p is close to 1, the PL 
lincshape also shows two maxima, however, there is a pe- 
culiar interplay in their magnitudes when the excitation 
energy is changed. This also matches the experimental 
data quite well. 

Such interplay can be easily understood. First (second) 
PL maximum corresponds to the distribution of ground- 
state energies of such subset of quantum dots, for which 
the lowest (highest) split-off level coincides with Ei aser . 
When Ei aser is below E m , the amounts of first type of 
dots is larger than that of the second type. When Ei aser 
is larger than E m , the dots of the second type win. 

The interplay is absent in PLE, because for each fixed 
ground-state energy there is always equal amount of low- 
est and highest split-off levels. 

The magnitudes of o\ and cr can be obtained indepen- 
dently as follows. For p — 0.94, the maximum of D p (e; a) 
lies at e = 0.765c. Then, from the position of the PLE 
maximum, eri=14.5 meV/0. 765=19 meV. The spread of 
the ground states Co can be determined independently 
from the nonselective PL [curve 9 in Fig. 6(a)]. It gives 
ao=18.2 meV. 

If we know the shape of the wave function, we may 
find the values for the parameters <7o,i and p, using Eqs. 
(2,4). As a guess, we may try the wave functions for the 
cylindrical quantum well with infinite walls: 

t/> m (r) « c os(^-)J m (^p-), (12) 
h R 

where h and R are correspondingly the thickness and 
the radius of the quantum dot, and v m is the root of 
the Bessel function J m . The integrals of interest, J mn — 
J d 3 r V'mV'n) are equal to: J 00 = 2.098, Jn = 1.552, and 
Joi = 1.435. 

It is easy to find all parameters in this approximation. 
First, let us estimate the magnitude of the spread. To 
find (To one has to know the volume of the quantum dot. 
Taking it to be the volume of the cylinder with the thick- 
ness 2.5 nm and diameter 25 nmp and using the param- 
eters of InGaAs: xo = 0.5, lattice constant a^p 0.585 nm 
and (dE g /dx) — 1.16 eV (both at x = 0.5) p we obtain 
(To = 13 meV. This value is less than the experimental 
value of 18.2 meV. It, however, shows that at least a sig- 
nificant part of the spread is caused by the composition 
fluctuations. There are, of course, some other reasons 
for spreading. The total spread, however, cannot be less 
than the calculated value. 

Two remaining dimcnsionless parameters are: p = 
0.795, and (ao/°~i) ~ 1.16. Though the ratio <Jo/<J\ is in 
a reasonable agreement with the experiment, the value of 
the correlation coefficient p is significantly less than the 
experimental value of 0.94. Note that the dimensionless 



parameters depend only on the form of the wave func- 
tions. Though it is possible to relate the discrepancy in 
p to the unknown shape of the real wave function, there 
is a more serious reason for this coefficient to be closer 
to 1. 

Among the other causes of spreading of energy levels, 
which are not taken into account in the presented the- 
ory, there are fluctuations of the shape of quantum dots. 
The distortion of the shape of a quantum dot, even when 
small, cannot be represented as a potential perturbation 
in the Schrodinger equation. It's effect on the positions 
and splitting of the energy level E m can, however, be 
described by an effective secular matrix SH of the same 
form as Eq. (6). The only difference is that the matrix 
elements u and x + iy are not described by the Eqs. (7) 
any more. Instead, they are determined by the integrals 
of the derivative of the unperturbed wave function at the 
boundary of the quantum dot .El 

Thus, the effect of the shape fluctuations is only in 
renormalizing the parameters o~o, o~ m , and p. It is im- 
portant, however, that for the case of pure shape fluctu- 
ations, the parameter p, defined as a correlation coeffi- 
cient between u and eo, is exactly equal to 1. The reason 
is that in an axially-symmetic quantum dot, the normal 
derivative of the wave function at the boundary is just 
a number rather than a function of coordinate. There- 
fore, there should be an exact relation between u and eo 
for each particular shape distortion. It is natural to as- 
sume that the effective value of p is closer to 1 when both 
mechanisms are involved, than it is in the case of pure 
composition fluctuations. 

VI. CONCLUSION 

In the present paper a simple theory is developed, 
which allows to describe selective photoluminescence 
data from an array of quantum dots with random pa- 
rameters. The theory explains large apparent Stokes-like 
shift between emission and absorption energies as the av- 
erage distance between the ground and first excited en- 
ergy levels in the dots. 

It is shown that existence of a random spread of the 
interlevel distances in the dots causes deviation of the 
positions of the maxima of the peaks seen in the PL and 
PLE spectra. Such deviation can make it difficult to 
determine the properties of the statistical distribution of 
energy levels in the array. It is suggested how the proper 
parameters of the statistical distribution may be obtained 
from the experimental data. 

The random shifts and splittings of energy levels 
caused by a white-noise random potential in the dots 
are studied. The density of states and mutual level-level 
distribution function are obtained for the case of axially 
symmetric quantum dot. The energy- level distribution 
and the resulting PLE and selectively-excited PL spec- 
tra appear to be close to that observed in the experiment 
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(see Fig. 6). It is shown that the major part of the spread 
of energies observed in the experiment originates from 
the random potential caused by the composition fluctu- 
ations. It is also suggested that the random fluctuations 
of the shape of the dots also contribute to the spread. 
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FIG. 1. A schematic diagram of the PLE and selec- 
tively-excited PL intensities when interlevel distances in all 
quantum dots are the same. Each curve is the spectrum 
taken at fixed excitation energy (PL) or fixed detection energy 
(PLE), which are denoted by pluses. 

FIG. 2. A schematic shape of the intensity as a function of 
the excitation and detection energies for the system with two 
excited states. The left "ridge" corresponds to the absorption 
by the first excited states of the dots. The "PL" and "PLE 
lines" in the inset are the loci of the positions of the maxima 
of intensity in the PL and PLE crossections respectively. The 
dashed lines show the lines of equal intensity. The maxima in 
a crossection occur at such points where the crossection lines 
(horizontal and vertical solid lines) are tangential to the lines 
of equal intensity. 

FIG. 3. The function Do(z;l) (full line), as given by Eq. 
(9), determines the density of states in the vicinity of the ex- 
cited level. The approximation by a Gaussian with an effective 
dispersion (dashed line) is also plotted. 
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FIG. 4. The two-level mutual distribution function 
P(eo,ei) given by Eq. (10) is plotted for fixed eo = (the 
"PLE" curve) and for different fixed ei (the set of "PL" 
curves). The curves are shifted upward arbitrarily to ease the 
reading. The correlation coefficient is: (a) p = 0.94 and (b) 
p=l. 

FIG. 5. The positions of the PL (crosses) and PLE (pluses) 
maxima taken from the spectra in Ref. 1 are shown in the 
combined plot. The resulting crossing lines make clear how 
the positions of maxima deviate from each other and from the 
major axis of the first "ridge" (dashed line). The crossection 
point is the maximum of the ridge. The second set of the 
PLE maxima is not shown because the corresponding peaks 
are seen not very clear in the spectra. A constant value of 
1290 meV has been subtracted from all energies. The dotted 
fine shows the unit slope, Ei aser = E detector . 

FIG. 6. (a) The experimental data from Ref. 1. replot- 
ted for comparison. The upper curve is the PLE spectrum 
obtained with the detection energy marked by an arrow. 
Pluses above the PLE spectrum show the fixed excitation en- 
ergies for each of the PL spectra. The curve 9 corresponds 
to the above-barrier excitation energy, (b) The PLE and se- 
lectively-excited PL lineshapes obtained from Eq. (10) with 
parameters p = 0.94, <ro = (J\ = 19 meV. Average energy 
positions of the ground and first excited energy levels are 
E = 1276 meV, and E\ = 1351 meV. 
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